MZ-TH/06-01 
PITHA 06/01 

February 1, 2006 



SUSY Ward identities for multi-gluon helicity 
amplitudes with massive quarks 

Christian Schwinn* 

Institut fur Theoretische Physik E 
RWTH Aachen, D - 52056 Aachen, Germany 

Stefan WeinzierF' 

Institut fur Physik, Johannes- Gutenberg- Universitat 
Staudmgerweg 7, D-55099 Mainz, Germany 



Abstract 
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involving a pair of massive quarks to amplitudes with massive scalars. This allows 
to use the recent results for scalar amplitudes with an arbitrary number of gluons 
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1 Introduction 



Recently, a number of efficient methods for the calculation of helicity amplitudes in QCD 
have been introduced, motivated by the relationship of QCD amplitudes to twistor string 
theory [1]. In the Cachazo - Svrcek - Witten (CSW) construction [2], tree level QCD 
amplitudes are constructed from vertices that are off-shell continuations of maximal helicity 
violating (MHV) amplitudes [3] , connected by scalar propagators. This formalism has been 
extended to massless quarks [4] and loop diagrams [5]. Subsequently a set of recursion 
relations has been found [6] that involve only on-shell amplitudes with shifted, complex 
external momenta. These methods extend earlier recursive formulations [7-9] which are 
based on currents with one off-shell leg. Recently a derivation of the CSW rules from the 
on-shell recursion relations has been given [10] and proofs of the on-shell recursion relations 
using conventional Feynman diagrams have appeared [11,12]. Another elegant technique 
applicable to scattering amplitudes with massless particles is the use of supersymmetric 
Ward identities (SWIs) [13-15] to relate helicity amplitudes involving massless quarks to 
those involving massless scalars or to purely gluonic ones [14,16-18]. 

In the light of the major role played by top quark physics at the LHC, it would be 
desirable to extend such methods also to multi-parton processes involving massive quarks. 
For instance, top-pair production with one or two additional jets contributes to the back- 
ground for Higgs production by vector boson fusion or in association with a top quark 
pair [19]. Tree amplitudes with a top pair and additional jets are also required for the 
next-to leading order corrections to top pair plus jet production [20]. Helicity amplitudes 
for a top pair plus three gluons can be found in [21]. Tree level helicity amplitudes with 
massive quarks are also needed within the unitarity method [22] for one-loop amplitudes, 
where a heavy quark is circulating in the loop. In this context, four point amplitudes have 
been computed in [23] and four and five point amplitudes have been obtained from on-shell 
recursion relations in [24]. 

The on-shell recursion relations are also applicable to massive particles [25, 26] and 
closed expressions for some sets of amplitudes with massive scalars and an arbitrary number 
of gluons have recently been found [27]. In the MHV approach it has been possible to 
include single external massive gauge bosons or Higgs bosons [28]. For massive quarks, 
the MHV approach cannot be used in a straightforward way since there are nonvanishing 
helicity amplitudes that go to zero in the massless limit where less than two partons have a 
helicity opposite to the remaining ones. In the on-shell approach only four point functions 
with massive quarks have been computed so far [26]. 

As a first step to extend the formalism of [2] towards massive quarks, we have reformu- 
lated in [29, 30] the Feynman rules of QCD including massive quarks so that only scalar 
propagators and a set of cubic and quartic primitive vertices without Lorentz or spinor 
indices appear. The cubic primitive vertices in this approach can also be used directly in 
the on-shell recursion relations. 

Recently, a particular set of multi-gluon helicity amplitudes with a pair of massive 
quarks has been calculated in [31]. These results should be related through supersym- 
metric relations to the one obtained in [27] for multi-gluon helicity amplitudes with a 
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pair of massive scalars. In this paper we work out the explicit transformation laws. The 
derivation follows the lines of the original paper [13]. There are several motivations for 
reconsidering the supersymmetric transformation laws: First of all, by taking into account 
recent developments in the definition of off-shell polarization spinors, we obtain short and 
elegant results. Secondly, the results allow us to relate amplitudes with massive scalars to 
amplitudes with massive quarks. Since there are more results available for all-multiplicity 
amplitudes with massive scalars, this is useful for top quark physics. Finally, we are very 
careful to keep internally all signs correct. This is an issue since there are different sign 
conventions in the supersymmetric and QCD community. The final result can of course 
easily be related to any other sign convention. Furthermore we have to respect the relative 
sign between amplitudes with anti-fermions and without. We verify the correctness of the 
transformation laws by a direct computation of amplitudes with a massive pair of quarks or 
scalars by using Berends-Giele type recursion relations and the diagrammatic rules of [29] . 

This article is organized as follows. In section 121 we review the notation for colour 
ordered partial amplitudes and the conventions for spinors and polarization vectors. In 
section El we introduce our supersymmetric model, give the supersymmetric transformation 
laws and apply them to derive relations among various amplitudes. In section0]we consider 
as an example amplitudes with a massive pair of quarks or scalars plus an arbitrary number 
of positive helicity gluons. Our conclusions are contained in section E3 

In appendix El we list our sign conventions. Appendix El contains the derivation of the 
supersymmetric transformation laws. In appendix [0 we extend the diagrammatic rules 
obtained in [29] to the supersymmetric model. In appendix [D] we give some technical 
details on the solution of the recurrence relation discussed in section HI 

2 Review of notation and off-shell continuation 

2.1 Colour decomposition 

In this work we are concerned with the evaluation of helicity amplitudes with one massive 
quark pair. We employ the usual decomposition of the full amplitude A n into gauge 
invariant partial amplitudes A n defined by [7,32,33]: 

Q n) 

g n-2 ^ (W)...TM«-i))..A n {Ql, g< 2),...,9,(n-l),Qn) , (2-1) 

creS„-2(2,...,n-l) 

where the sum is over all permutations of the external gluon legs. The colour structure 
is contained in group theoretical factors given by traces over the generators T a of the 
fundamental representation of SU(3). 

We will also encounter amplitudes with scalars transforming under the fundamental 
representation of SU(3). The colour decomposition is similar to the one in eq. (j2.1|) : 

An((f>i, g 2 , g 3 , ...,g„-i,(p n ) = 
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g n-2 J- (T a ^...T a ^%A n {h,9«(2),..,g« { n-l)An) ■ (2.2) 

<r€S„_ 2 (2,...,n-l) 

2.2 Polarization vectors and spinors 

It is a well known fact that any four- vector can be decomposed into a sum of two light-like 
four-vectors. Given a fixed light-like four-vector q, one can associate to any four-vector k 
(which need not be light-like) a light-like four- vector k , defined by 

This prescription has also been used for the off-shell continuation of MHV amplitudes [34] 
in the context of the CSW rules. Unless stated otherwise, we will use the same reference 
vector q for all momenta. The projected four-vector k b can be used to define continuations 
off the null-space of spinor products as follows: 

(hh) = (k[ - \k 2 +) , \k 1 k 2 ] = {k\ + \k\-). (2.4) 

The two two-component spinors \k b +) and \k b — ) are the Weyl spinors corresponding to 
the projected momentum k bfl . These can be used to define the off-shell continuation of the 
polarization vectors for gluons: 

# M) = ± Mla!M . (2 .5) 

Helicity methods can be extended to massive fermions [17,31,35,36]. We introduce an off- 
shell continuation of massive spinors, using the projection ()2.3j) as for the gluon polarization 
vectors: 



ft + m 

u(k,±) = (q T \ 



(qT\k b ±) 



v(k,±) = (qT\ ' 





m 


l?±>, 


(k b =F |?±) 


(k b ±\ + 


m 


(1T\, 


(qT\k b ±) 


l* b =F> " 


m 


\q±), 


(k b =F |?±) 


(A; b ±| - 


m 




(qT\k b ±) 



(2.6) 



The normalization is chosen in order to allow for a smooth massless limit. We label 
the helicities as if all particles were outgoing. As a consequence, the spinors u{k) and 
v(k), which correspond to particles with incoming momentum, have a reversed helicity 
assignment. 

The reference momentum used in the definition ()2.6|) is not an unphysical quantity 
that has to drop out in the final result for the helicity amplitudes, as in the case of the 
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definition of the gluon polarization vectors, but rather defines the quantization axis of the 
quark spin [35]. Since for any given choice of q the spinors ()2.6)1 form a complete basis of 
solutions of the Dirac equation, the calculation can be simplified without loss of generality 
by a suitable choice of q. The amplitudes for a desired physical polarization can be obtained 
by a straightforward linear transformation. 



3 SUSY Ward identities for massive quarks 

Calculating examples of helicity amplitudes for a massive quark pair and positive helicity 
gluons, one observes a simple relation to amplitudes with massive scalars that have been 
calculated in [25,27,37]. As a simple example, we find the four point function 

A(Q+ a + a+ 0-) - {Aq) [23] W (3 1) 

The results for an arbitrary number of positive helicity gluons are quoted in appendix [01 
Comparing to the scattering amplitude for two massive scalars and two positive helicity 
gluons [25,27,37] 

4/7+ + i \ [23] 2im 2 . , 

A <*-» + -**-*> = M (*. +*.)'-*,■ ■ (3) 

we see these amplitudes are related by a factor similar to the SUSY relation between 
quark and gluon amplitudes. In the next subsection we define a super symmetric toy model, 
which contains QCD as a subset, and derive the supersymmetric transformation laws. In 
subsection 13.31 we apply these relations to amplitudes involving a pair of massive quarks 
or scalars. It turns out, that the relations are particular simple, if the reference spinor 
defining the spin axis of the heavy quarks, is chosen to be same for both the quark and the 
anti-quark. In subsection 13.41 we discuss modifications which arise if the reference spinors 
are not identical. 



3.1 Definition of the model 



To embed QCD with massive quarks into a supersymmetric extension, we have to use two 
chiral super-multiplets $ + = (<p + ,ip+,F + ) and $_ = (<^_,^_,F_) and combine the spinor 
component fields into Dirac spinors as follows 

% =(*/>-,$+), *=($+)■ ( 3 - 3 ) 
The chiral fields are written in component fields as follows: 

$ + {y, 9) = V+ (y) + y/2fol> + {y) + 9 2 F + (y), y» = - i9a"9, 

0) = <p_(y) + V29ip-(y) + 9 2 F„(y), y» = x» - i9a<*9. (3.4) 
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The Lagrange density of our model is given by 

C = ^-Tr WW\ F + ^-Tr WW\ p 
8g 2 8g 2 lF 



D 



D 



+ m $_$ + | F + m $ + $_| (3.5) 



where the vector multiplet V a = (A a , X a , D a ) contains the gluon A a and the gluino A a . 
The field strength is defined by 

W A = -~ {DD) (e 29V D A e- 29V ) , W A = ~ ( DD ) ( {D A e ~ 29V ) e 2gV ) , V = T a V\ (3.6) 

and V a is given in the Wess-Zumino gauge by 

V a = (9^9) Al + i9 2 (9\ a ) - i9 2 (9\ a ) + ^9 2 9 2 D a . (3.7) 

For the gluino field we define a four-component spinor as follows: 



A = -i\ A ) , A = 7^ . (3. 



It is also convenient to redefine the scalar fields as follows: 

<P- = <f-, 0+ = <P-, 0+ = <P+, 0- = <P+- (3-9) 
Eliminating the D- and F-terms and with the convention that 

$±y = 4> T (3-io) 

we obtain 

C = -- (F^) 2 + % -h a ^Dfh b + zf rfDJi! - (3.11) 

+ (D M + ) t (Dy + ) - m 2 0„0 + + (D^) j (D"0_) - m 2 + 0_ 
-V2g [ ( j) + A a T a P + m + ^P_A a T>_ - 0_A a T a P_^ - ^P+A a T a 0+] 

(0 + T a 0_-0_T a + ) 2 , 

where the covariant derivative in the fundamental and adjoint representation is given by 

D, = d, - igT a A*, Df = d, - gf^Al (3.12) 

The chiral projectors are denoted as usual by P± = 1/2(1 ± 75). As in [29] we can derive 
primitive vertices for the additional interactions that arise compared to nonsupersymmetric 
QCD. The results are collected in appendix IU1 
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3.2 Transformation of the fields 



The SUSY transformations of helicity states have been first derived in [13] and later applied 
to QCD helicity amplitudes [14,16]. We use a commuting SUSY generator Qsusy(?7) = 
V A Qa + VaQ A obtained by multiplying the SUSY generators by a Grassmann valued spinor 
rj. This generator satisfies the SUSY algebra 

[QsusyM, Qsusy(0] = ( 3 - 13 ) 

SUSY transformations of generic bosonic or fermionic fields $ are generated by 5 V ^ = 
[Qsusy(^), ^]- With the notation 

rJ(*) = V2fo±|*=F>, ^(k,q) = V2m^0 ) (3.14) 

one finds the transformation laws for outgoing particles: 

5 V Q± = -rJ(A;)0 ± + E*(k, qW, 5^ = -V*{k)Q± - Z*(k, q)Q*. (3.15) 

Details on the derivation are given in appendix [Bj Eq. ()3.15|) is the main result of this 
paper and generalizes the well-known results for massless particles to the massive case. 
Compared to the massless case there are additional terms with T,*, which are proportional 
to the mass of the particles in the super- multiplet. For practical applications, one can 
achieve a considerable simplification by noting that q) = 0. Therefore the additional 
contributions to ()3.15|) can be avoided provided the SUSY spinors |r/±) are taken to be 
equal to a Grassmann parameter 9 times the reference spinors |g±) used in the definition 
of the massive spinors (j2.6j) . This requires that the same reference spinors |g±) are chosen 
for all massive quarks in the amplitude. In this case we have the following transformation 
laws, which are identical with the massless case: 

8 n Q± = -T±(k)4>±, 5^ = -T*{k)Q±. (3.16) 
The gluon and the gluino transform as follows: 

8 v g± = -rJWA* = rJ(A;)A ± , 

^A ± = -rj(% ± , ^A ± = r^(% ± . (3.17) 

3.3 SUSY Ward identities for massive quarks and scalars 

We can now derive relations connecting amplitudes of massive quarks to amplitudes of 
massive scalars, generalizing the example of the four point amplitudes (J3.1|) and (|3.2|) . 
In this section we will always choose rj = 9q and use the same reference vector q in the 
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definition of the spinors for all massive quarks. Therefore terms proportional to S ± do not 
contribute. Furthermore we treat in this section all particles as outgoing. 

Since the SUSY charge Qsusy annihilates the vacuum for unbroken SUSY, we obtain 
SUSY Ward Identities (SWIs) from the relation [13] 

= (0|[Q S usy(*7), zi...z n )\0) = ^2(0\z 1 ... [Qsvsriv), 4 ■ ■ ■ z n \0) , (3.18) 

i 

where the Zi are arbitrary creation or annihilation operators. 
Applying the SUSY transformation to the amplitude 

A(^t,gt,-,9n-l,Qn) (3-19) 

with one + scalar, (n — 2) positive helicity gluons and a negative helicity quark results in 
the relation 

n-l 

r-ikMiQf, 9t, gtv Qn) ~ E r v( k ^t, at, ■ ■ ■ , A+, • • • , gtv Qn) 



i=2 



r-{k n )A(4>+, gt, ■ gti, <f>n) = o. (3.20) 



Note that these identities hold to all orders in perturbation theory in the supersymmetric 
toy model. However, only on tree level the scalar and quark amplitudes are identical to 
those in pure QCD while in higher orders loops involving SUSY particles contribute. 

The terms arising from the transformations of the gluons vanish at tree level since there 
are no primitive vertices of the form V(<f)f,Af,Q~). Therefore we obtain the relation 

A n (Q+, g+, . . . , g+^, Q~) = £jU»(#, 9% \ • • • , gt-v C), (3-21) 

which generalizes the relation among the four point amplitudes ()3.1j) and (j3.2j) to an ar- 
bitrary number of positive helicity gluons. In the same way, the SUSY transformation of 
>g£ ) •••) <?n-i; ^t, ) gives the corresponding relation for exchanged quark helicities: 

A n (Q^gt, ...,g+_ lt Q+) = -M^fc <? 2 + , g+_ 1: 0+). (3.22) 

We can also derive SWIs for amplitudes including a negative helicity gluon. Leaving 
the helicity of one massive quark arbitrary, these relations read 

o = 8 v A{4>t, gt, ■■■■u J Q$) r-fa^Qf, gt, ■ ■ ■ ,gj,. ■ ■ ,Qt) 

- r-(kj)A$t, !h \ Qn) - r±(k n )A$+, !h g] 4>t), (3.23a) 

where the terms from the transformation of the positive helicity gluons vanish for the same 
reason as above. For the amplitude with two positive helicity quarks A(Qf, . . . , g~ , . . . , Qn), 
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the last term vanishes since it involves two different scalars <p + and (f) + . This 'helicity flip' 
amplitude is therefore directly related to an amplitude involving a gluino and a scalar. 

A(Ql,gt,...,gJ,...,Q+) = ^A(ft,g+, . . . , Aj, . . . , Q+). (3.23b) 

For the amplitude with one negative helicity quark, all three contributions survive. The 
corresponding identity for the opposite quark helicities is given by 

= S v A(Qf, g+, ...,g~,..., 0+) = F^(ki)A{4t, gt, .. .,gj,. . .,</>+) 

- T-(k 3 )A(Qf, g+,...,Aj,..., <f>+) + F-(k n )A(Qf , g+, . . . ,gj, . . . , Q+). (3.23c) 

We can simplify the identities (|3.23|) further by a specific choice of the reference spinors. 
With the choice \q+) = \j+) we obtain 

, , . , = 0, (3.24) 

\q+)=\3+) ' V ' 

\ q+) =\ j+ ) = J[j} An ^> 9t,---,9j,---, <Pn), (3-25) 

k+H ;+> = ~^j) An ^ T ' ° ;/h (3 ' 26) 

The scalar amplitudes appearing in these identities have been computed for an arbitrary 
number of gluons in [27] using on-shell recursion relations. However it should be noted 
that the choice \q+) = \j+) defines the quantization axis of the spins of the heavy quarks. 

Finally we can consider amplitudes with two negative helicity quarks and positive he- 
licity gluons only. In this case the SWI relates the quark amplitude to a sum of gluino 
amplitudes: 

= 5 n A(^,g+, ...,g+_ v Q~) = T^kMiQi , at , 9tv Qn) 

n-l 

- r+(^)A(0r, s 2 + , V <?+_!, Qn), (3.27) 



A(Qt 


,gt> ■ 


■ ■,gj,- 




A(Qt 


,gt, ■ 




■ ■ ) Qn 


A(Qi 


,gt, ■ 


■ ■,gj,- 


■ ■ ) Qn 



i=2 



where the additional term from the transformation of the quark vanishes. 

These last examples show that SUSY methods can be extended to amplitudes with 
additional negative helicity partons. However in general these relations will involve ampli- 
tudes with external gluinos. 

3.4 Alternative choice of reference spinors 

In [31] a different choice of auxiliary spinors than in (|2.fij) has been used to calculate helicity 
amplitudes for a massive quark pair and an arbitrary number of positive helicity gluons 
using the Berends-Giele relations. In this subsection we show how to relate our results to 
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those of [31]. In that work the momenta of the massive quarks k\ and k n are decomposed 
in terms of two lightlike vectors ky n according to 

fc? = (3 + k» + (3_K, K = p_h* + p + %£, (3.28) 



with f3± = |(1 ± (3) where (3 = J 1 — p is the quark velocity. The spinors are 

then defined similarly to ()2.6|) but using (3 + k n as reference momentum for Qi and (3 + k\ as 
reference momentum for Q n . Explicitly they read: 



fl(1 - ±> = <nTl = "+ (1±l + WW^) (nTl 



(3.29) 

in ^ 

Pl /2 (n T |T±) " + ' ' $T (n =F |T±> ' 



Since the choice of different reference spinors for the two quarks corresponds to different 
spin-vectors, different helicity combinations give a nonvanishing result compared to our 
formalism where the same reference spinor is used for both quarks. The general trans- 
formations ()3.15j) still relate amplitudes with massive quarks to amplitudes with massive 
scalars. As an example we consider the SWI obtained from 

6r,A$t,9t,~,9t-i,Q$) = 0. (3.30) 
The explicit forms of the relevant transformation formulae from ()3.15j) read: 

T- |r/+ X 



\ I (l 



. 71' 

\n+ 



(n - \r)+ ) _ 



5 n <t>i = V2/3^(t]- \i+)Qt + V2mp- 1/2 ) — ^4or> (3-3i; 

x 1 (n-\l+) 

and by choosing r] = 6k n we obtain for the helicity flip amplitude 

— TTl — 

A(Q+,g+,...,g+_ l ,Qt) = tz ^A(<j>+,g+,...,g+_ 1 , ( f>-). (3.32) 

/3+ (l - |n+) 

This relation can immediately be verified for n = 4: Ref [31] obtains for the four point 
function the result 

A (Qi, gt-> dt i Qt ) — zz ~ \ tt^Jtt; , r \ 2 II?- (3.33) 

(3 + (l - |4+) ( 23 ) (*i + k v m 

The scalar amplitude is given by [27] 

Hit, 9 i,9tAl)= [ § ) W 0^ . (3.34) 
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4 Recursion relation and results for amplitudes with positive 
helicity gluons 

In this section we consider the degree zero amplitudes 

A(Q+,gt,...,gt-i,Qn) and M$t,9t, ■ ■ ■ ,9n-i, fa) ( 4 - x ) 

with a massive pair of quarks or scalars plus an arbitrary number of positive helicity gluons. 
In the following we will discuss amplitudes with one particle off-shell, this particle will be 
marked by a hat. It is also convenient to treat Q n as an incoming quark rather than an 
outgoing anti-quark. Let us take the outgoing particles 1 to n — 1 with positive energy. 
Therefore p n = k\ + ... + fc n _i = —k n has a positive energy component. With our phase 
conventions, the corresponding amplitudes are related by 

A(Q+,g+,...,g+_ v Q^) = iA(Q+, g£,..., £+_il<2p n =-fcJ- (4-2) 

In the following we will use the diagrammatic rules of [29] that are briefly reviewed in 
appendix O The amplitudes in eq. (j4.1j) are of degree zero, therefore in our conventions no 
helicity-flip vertex occurs. As a consequence, the corresponding amplitudes with one-leg 
off-shell satisfy rather simple Berends-Giele type recurrence relations: 

n-l 

j=2 1,3 

-^-A n _ j+1 (gf,...,g+_ l ,g_ {jtn) ), (4.3) 

where we have defined 

3-1 

l=i 

The corresponding recurrence relation for the scalars reads: 

n-l 

A n (fa, g+_ x , J") = V 3 (<l>tj, 9j, w fa) ij2— — a ->9j-v fan j)) 

j=2 ij 

A n _ j+1 (g+, ...,g+_ x ,g_^ n) ). (4.5) 

^> 

Here the quartic primitive scalar-gluon vertices do not contribute since the off-shell con- 
tinuation of the polarization vectors (|2.5|) implies e + (k iy q) ■ e + (kj, q) = 0. 

The recursion starts with the two-point amplitudes, which are given by the inverse 
propagators: 

A 2 (g+,gZ j ) = -ikl A 2 (Qt\Q 1 ) = A 2 (faXi) =-i{k\-m 2 ). (4.6) 
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These recurrence relations are solved explicitly in appendix [Dj In the present example we 
can use the recursion relation to establish the relation ()3.21)1 between the amplitudes of 
quarks and scalars directly. 

After applying crossing symmetry the quark amplitude entering the SWI can be written 

as 

A n (Qt, ■■■i9n-iiQn) = iAnlQi, gt-\\Q Pn =-k n ) 

n-l 



i 

T 

j.n 



— ^ % 

-«) A i (Qt, ■ ■ ■ , 9j- 1 , QZ (u) ) p- At-j+i (gf , • • • , 9n- 1 , 9-(j, n ) ) • ( 4 - 7 ) 



We now use the induction assumption that the identity ()3.21j) holds up to n — 1 already 
for the one- leg off-shell amplitudes. This is true for n = 2: 

A 2 (Qf,Q 2 ) = iASt,fc)- (4-8) 

Since the primitive cubic scalar-gluon vertex (jC.3|) and the quark gluon vertex (jC.l|) satisfy 
the identity 

v(Qt, i + , Q; n= -j = ^v$t, i\ c) = -^n$i, i + , <pz), (4.9) 

we find that every term of the sum in the recursion relation contains a product of the form 

V 3 (Qtj,9tn, QZ n )MQ+, QZ {1>j) ) 



Tj^vStj, gin, <t>n)M4>t, ■ ■ ■ U - fcd,)) (4.io) 

W) 

sum. The remaining factors give the scalar amplitude (|4.5jl and we obtain the SWI (j3.21j) : 



Therefore in every term appears the same overall-factor that can be pulled out of the 



A(Qt,gt, ...,9ti,Qn) = j^A($t,9t,...,9ti,fa)- (4.11) 

This is a non-trivial cross-check on the correctness of the supersymmetric transformation 
laws in eq. (j3.15|) . in particular since the equivalence of the analytic expressions for the 
all-multiplicity quark amplitudes [31] with the scalar amplitudes [27] is not immediate. 
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5 Conclusions 



In this paper we considered relations based on supersymmetry between multi-gluon am- 
plitudes with a pair massive quarks and corresponding amplitudes where the quarks are 
replaced by massive scalars. These relations generalize the well-known formulae for the 
massless case. The relations are useful as they allow to obtain amplitudes with heavy 
quarks from amplitudes with scalars. Furthermore, if both sets of amplitudes are known, 
the super symmetric relations provide a valuable cross-check. 

Note added in proof 

After the submission of this paper, a form of the scalar and quark amplitudes with an 
arbitrary number of positive helicity gluons has been found [40] that is much simpler 
compared to those given in [27, 31] and in (JD.7)) and that has been used to verify the 
identity 1)3.32)1 explicitly for an arbitrary number of gluons. 
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A Notation and conventions 

In this appendix we give a comprehensive summary of our notation and sign conventions. 
The convention for the metric tensor is g^ v = diag(+l, — 1, — 1, — 1). We use the Weyl 
representation for the Dirac matrices 



where a^ B — (1, —a), a^ AB = (1, a) and o = (cr x , a y) a z ) are the Pauli matrices. The chiral 
projectors are as usual P± = 1/2(1 ± 75). The two-dimensional antisymmetric tensor is 
defined by 




(A.l) 




(A.2) 



Indices of two-component Weyl spinors are raised and lowered as follows: 



= gAB^ ^A = £ AB^^ ^ = ^A £ ^ b = ^A £ab ^ 



(A.3) 



The generators of supersymmetry are taken to be 



(A.4) 



'.4 
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They satisfy the supersymmetry algebra: 

{Qa,Q A } = -2< A ^- (A.5) 

Spinor products are defined as 

= Z A VA, £f} = £ A fj A . (A.6) 
In the bra-ket notation spinor products are denoted as 

(pq) = (p- \q+) = p A q A , [qp] = (q + \p-) = q A p A , (A.7) 
In terms of the light-cone components 

P+=P0+P3, P-=P0~P3, P±=Pl+iP2, P±*=Pl~iP2- (A.8) 

of a null- vector p^, the corresponding massless spinors (p ± | and |p±) can be chosen as 

e l 2 f - Pj _, \ e ! 2 / p^ 



-it -i± 
e 2 e 2 

{P+\ = —f^=r{-P±,P+) , {p-\ = —n=^ (P+'P±V ' ( A - 9 ) 

vb+l vb+l 

where the phase (f) is given by 

= |p+|e*. (A.10) 

If p + is real and p + > we have the following relations between a spinor corresponding to 
a vector p and a spinor corresponding to a vector (— p): 

|(-p)±> = i\p±), 

((-p)±\ = i{p±\. (A.ll) 

Therefore the spinors of massive quarks and anti-quarks are related by u(—k, ±) = iv(k, ±) 
and u(—k, ±) = iv(k, ±). The polarization vectors of the gluons are unchanged under the 
reversal of the momentum. Denoting a partial amplitude with n incoming particles with 
momenta p, and m outgoing particles with momenta kj by 

A n+m ($ fcl ...$ fc ,J$ Pl ...$ p J (A.12) 

the amplitude with an outgoing anti-quark is related to that with an incoming quark by 

A{.., Q k]) ...\...) = -iA{...\..., Q Pl= . kj , ...) (A.13) 

where the notation used for the quarks follows that of the spinors, i.e. incoming quarks 
and outgoing anti-quarks are denoted by Q and outgoing quarks and incoming anti-quarks 
by Q. The same phase as in (jA.13|) appears for the exchange of an outgoing quark with 
an incoming anti-quark while for gluons and scalars no phase appears and it is sufficient 
to reverse the momentum. 
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B SUSY transformations of helicity states 

In this appendix we give the derivation of the supersymmetric transformations in eq. (j3.15j) 
and eq. (j3.17J) . We follow closely the approach of ref. [13]. The Lagrangian in eq. ()3.5|) 
is invariant under supersymmetric transformations. The supersymmetric transformations 
act on the component fields of $ + as follows: 

5 v ip + = V2rjif: + , 

5 v ip+ = -V2ia"fjd tl (p + + V2r]F + , 

5 V F + = -V2i^d^ + . (B.l) 
Here we used the two-component Weyl spinor notation. The components of $_ transform 



as 



8 v (p- = \/2#-, 

5 V F_ = -V^irja^d^-. (B.2) 
Finally, the components of the vector multiplet transform as 

8 v D a = rja^d^\ a + d fl \ a cr IJ 'fj. (B.3) 

We recall the connection between the two-component Weyl spinor notation and the four- 
component Dirac spinor notation: 



A = (iX A , -fX A ) , A = ( j . (B.4) 

For the scalars we used the redefinition 

4>- = V 5 -, 0+ = <P-, 0+ = <fi- = <P+- (B.5) 
together with the convention that 

(4>±y = <h- ( B - 6 ) 

The asymptotic fields we may expand in terms of creation and annihilation operators: For 
the scalars and the quarks we have: 
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Let us consider all particles in an amplitude as out-going. To derive supersymmetric 
relations among different amplitudes with out-going particles we have to know through the 
reduction formulae the transformation laws of the annihilation operators, e.g 

AiQHh), Q ± (h), ± (%), <f>±(k j+ i), ■■■) = 

(0\b{h, ±), d{k 2 , ±), . . . , a±(%), a±{k j+1 ), . . . |0> . (B.8) 



The annihilation operators can be projected out in eq. (|B.7|) : With the notation 

fd,g = f(d,g)-(d,f)g (B.9) 

we have 

o^Jfe) =i J d 3 x (e ikx d , a±(k)=i J d 3 x (e ikx d T (x)) , 

b(k,\) = [ d 3 xe ikx u{k,X)j°^{x), d{k,X)= [ d 3 xe ihx ^{x)^°v{k,X). (B.10) 



These equations allow us to obtain the transformation laws of the annihilation operators 
from the known transformation laws of the fields. For example 



5 v a_{k) = i J d 3 x [e ikx d 5„0+(x)) = ijd 3 x [e ikx d <^-(z)) 

= ,/*, (e- S VS*w) = i J d 3 x (e- S ^.»(,)) . (B.11) 
Reinserting the expansion eq. (jB.7|) for one finds the transformation laws 
5 v a ± (k) = ^d(M) (V^, -A)P± ( , 
<$„a±(A;) = ^(v / 2(r ? ,r ? )P ± n(A;,-A))6(A;,A). (B.12) 

A 

For the fermion field we find 

s A k,x) = -V&<M) ( ) . 

8r,d{k,X) = V2(a_(k)ri,a + (k)f])v(k,\). (B.13) 

Using the explicit quark spinors in eq. (|2.6|) and taking the sign change for amplitudes with 
anti-quarks discussed at the end of this section into account we obtain eq. (j3.15j) . Similar 
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considerations apply to the particles in the vector multiplet. The gluino is a Majorana 
fermion with expansion 

A(X) = ^ J (2^82fcO ( 6 ^> ~^ e ~ lkX + ^ A )^> A ) e ^) • ( B ' 14 ) 
The asymptotic expansion of the gluon reads 

A"(x) = J2J ^^(aC^^e^^^^ + atC^^e^^^e^). (B.15) 

The annihilation operator for the gluon is projected out by 

a{k,X) = -i [ d 3 xe ikx d £»{k,\)A>*{x). (B.16) 



We then obtain 

MM) = ^(k,X)[(r]-\r\k-)b{k,-)-(r] + \r\k+}b{k,+)) 



S v b{k,±) 



±u(k, X)YP± (j ) E A >(*, A ')' (B-17) 



There is one subtlety in writing down the transformation laws for the fields: For the 
application towards supersymmetric Ward identities we have to keep the relative sign 
between different amplitudes correct. Let us consider a Feynman graph corresponding to 
the matrix element 

(o \b...d...T [&riV>i...i?„r n v>J tf..J\ o) . (B.i8) 

Here ipjTjipj denotes a generic fermion interaction vertex. In order to apply the Feynman 
rules, we have to anti-commute the fields and the creation and annihilation operators such 
that 

(0 \ip b\k, A) | 0) -> u(k, A), (0 \b(k, A) $| 0) -> u(k, A), 
<0|^rf f (fc,A)|0) -> u(Jfe,A), (0|d(Jfe,A) V|0) ->u(M). ( B -!9) 

It is easily seen that this reordering brings a minus sign for each spinor v and w [38]. For 
an individual amplitude these signs are an overall factor and are usually ignored. However, 
if one relates amplitudes with different particle contents through supersymmetric Ward 
identities to each other, these signs have to be taken into account. We can incorporate these 
signs into the transformation laws, by adding an additional sign to each transformation 
law, for which the number of anti-fermions changes by one unit. For example 

<^± = -r±(A;)A ± , 8^ = T±(fc)A ± . (B.20) 

In the second equation we inserted an additional sign with respect to eq. (jB.17|) . which 
compensates for the change in the number of anti-fermions. 
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C Diagrammatic rules 



In [29] we have obtained diagrammatic rules involving only scalar propagators i/k 2 for 
gluons and i/(k 2 — m 2 ) for massive quarks. We introduce so called primitive vertices that 
are obtained by contracting the vertices of the standard colour ordered Feynman rules with 
the off-shell polarization vectors and spinors ()2.5|) and ([2.6)1 . Within this approach it is 
convenient to take the momentum flow of a fermion always in the direction of the fermion 
arrow line. That implies that we replace an outgoing anti-fermion by an incoming fermion. 
All other particles remain outgoing. The tri-valent primitive vertices of massive quarks 
include vertices present both for massless and massive quarks like 

V 3 (Ql,Q^gt) =<v^S V a {Q;,q£,to) = i>/2@jL. (C.l) 

In addition, there are primitive vertices involving a helicity flip along the quark line that 
vanish for massless quarks: 

V 3 (Qt, Qt, = iyftm J^j, V 3 (Qr, Qt, gf) = ~ l ^ m ^l^iy (C ' 2) 

One can define the degree of a vertex or of an amplitude as the number of "-"-labels minus 
one. In the diagrammatic rules, only primitive vertices of degree zero and one occur. 
Furthermore, the degree of an amplitude is exactly the sum of the degrees of the primitive 
vertices. 

It is also straightforward to obtain the primitive vertices involving the scalars and 
gluinos resulting from the Lagrangian (j^.lljl . For the scalar interactions with gluons we 
obtain the primitive tri-valent vertices 

V(tt,<t>l,9i)=iV2 [ ^, Vifafagj) =-^^1^. (C.3) 

The primitive vertices resulting from the couplings between quarks, gluinos and scalars are 
obtained by inserting the spinors of the massive quarks ([2.6)1 . One finds the nonvanishing 
tri-valent vertices 

V(Q+, fa, A 3 + ) = -iV2[13), V(Q-[, <f>t, A 3 ) = iV2 (13) , 

V(4>t,Q2,As) = -iV2(32), V(^,Q+,Af)=iV2[32], 
V$+,Q+,As) = -iV^m|||, y(0r,Q 2 ",A 3 + ) = tv^m^|. (C.4) 

In addition there are four-valent vertices, which however we do not list here. As for 
pure QCD only degree one and zero vertices occur. For the application in the SWIs it is 
important to note that the scalars only couple to gluinos with the opposite helicity label. 
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D Solution of the recursion relation 

In this appendix we give some details on the solution of the recurrence relations (J4.3)) 
and (|4.5|) . We follow closely the calculation of [31]. For the solution to the recursion 
relation (j4.3J) we make the ansatz 

Pn [(m + k V ~ m \ L J 

where A n (g+, . . . , gn-i\% n ) is the amplitude for massless quarks with one leg off-shell [7,32] 

Mit, • • - ^-'H^) (12> ... ((K _ 2 ^ 1))(( „_ 1)9> - (D.2) 

In the massless limit, the formula ()D.1|) reduces to the known result (jD.2|) while in the 
on-shell limit a finite term proportional to m 2 remains that is determined by the quantity 

A(Q+ 10- ) im22n '" 1 MB n 



The structure of the amplitude (jD.lj) is similar to the current with two off-shell gluons 
obtained in [39]. 

In the recursion relation (J4.3)) also enters the well known expression for the gluon 
amplitude with one negative helicity gluon [7,8,32,33] and one off-shell leg 

A n {gt,...,gU,9n) = Tl 2 ~\-ik\ n ) { ^ (12) - {{n _ 1}) {{n _ . (D.4) 

As a check, we have reproduced this amplitude from the Berends-Giele relations within 
the formalism of [29]. 

Inserting the ansatz (jD.l|) into the recursion relation ()4.3|) . one obtains after some 
tedious steps a recursion relation for B n that is equivalent to that considered in [31], 
although there the part vanishing on-shell has been split off differently. Adopting the 
solution given there, we further split B n into two contributions 

(12) 



B,, 



(lq) (fcf 4 - m 2 



{k I n -m I ){2+\Uq+)+B n , (D.5) 



with B A = 0. The function B„ satisfies the recursion relation 



'n 

n-l 



B n = (2 + \hh,M+) + £ ^P^H^TF^ *V < D ' 6 » 



A closed solution to this equation can be found in [31]. For n = 4 we therefore obtain the 
result quoted in (|3.1j) while the on-shell amplitudes for n > 5 are given by 



im 



2 n / 2 - 1 (nq) B n 



A n (Q+, ...,g+_i, Qn) (jfe 2 g _ m2)(A ,2 4 _ m2) (lg) (23) _ _ ^ n _ 1)g) ■ (D.7) 
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The explicit results for the five and six-point functions obtained from (jD.7|) and (jD.6|) are 
given by 

MQUUUtQl) - (D . 8 a) 

MQt,gt, 9 t,9t,9t,Q^ = 2V {lq) {5q) (23) (34>(45> ( fc 2 3 _ m2)(fc 2 4 _ m2) P.8b) 

(g-|Mfe|4+) (2+|^ 2)4 |4-) 



(2 + | ^2,^4,6 1 g+) - 



(fc? 5 - m 2 



Comparing to the results for scalars obtained in [25,27,37] one finds that the five point 
functions obviously satisfy the SWI (j3.21j) . while the agreement of the six point functions 
can be established after some use of momentum conservation and Dirac algebra. 
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